Structure and three-body decay of 9 Be resonances 
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The complex-rotated hyperspherical adiabatic method is used to study the decay of low-lying 
9 Be resonances into one neutron and two a-particles. We investigate the six resonances above the 
break-up threshold and below 6 MeV: 1/2 , 3/2 ± and 5/2 ± . The short-distance properties of each 
resonance are studied, and the different angular momentum and parity configurations of the 8 Be 
and 5 He two-body substructures are determined. We compute the branching ratio for sequential 
decay via the 8 Be ground state which qualitatively is consistent with measurements. We extract the 
momentum distributions after decay directly into the three-body continuum from the large-distance 
asymptotic structures. The kinematically complete results are presented as Dalitz plots as well as 
projections on given neutron and a-energy. The distributions are discussed and in most cases found 
to agree with available experimental data. 



PACS numbers: 21.45.-v, 21.60.Jz, 25.70.Ef, 27.20.+n 
I. INTRODUCTION 

The structure of 9 Be has been extensively studied, 
both theoretically [l|-[T6| and experimentally (l7r{27j , but 
there are still large uncertainties in the structure and de- 
cay of the low-lying excited states. This is surprising 
and worrisome in view of the large efforts and the ex- 
pected rather accurate approximation as a simple three- 
body system where the intrinsic degrees of freedom are 
inactive. Are the problems related to inaccuracies of the 
theoretical models, the numerical techniques, direct ex- 
perimental uncertainties, data analysis, or interpretation 
of the data in comparison with model results? 

In theory the three-body continuum problem is better 
handled and more accurately solved for nuclear systems 
with the special mixture of short and long-range interac- 
tions. Observables rather close to the directly measured 
quantities can be delivered. In experiments both beam 
quality, detector systems and systematic analyses have 
improved substantially in recent years. This means that 
genuine three-body systems can be treated fully and pre- 
cisely in both theory and experiment, more specifically 
complete kinematics of the fragments are available. The 
road to detailed comparison is therefore paved. The sim- 
plest systems should then be understood before reliability 
can be expected for more complicated scenarios. 

Furthermore, the results of requested applications in 
astrophysics, where often the energies are too low to be 
reached experimentally, can only be indirectly tested by 
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their implications. The approximations employed so far 
in predictions should then be tested by comparison. A 
reasonable procedure is to select a three-body system, 
compute and measure the best we can, and compare as 
detailed as possible. The choice of 9 Be is tempting as a 
rather simple system which is accessible to both theory 
and experiments. In addition, this is a system of partic- 
ular interest in astrophysics, where formation of 9 Be can 
proceed through the reaction a(cm, 7) 9 Be. The subse- 
quent reactions, 9 Be(a,n) 12 C, link to heavier elements in 
stellar nuclear synthesis responsible for the present Uni- 
verse. 

From the early days of Nuclear Physics, the structure 
of the 9 Be nucleus has been considered a prototype of the 
cluster-like structure of nuclei. Therefore, different types 
of three-body descriptions have been used to describe it: 
early cluster models [D-d, 0, [|[ , and more sophisticated 
ones, e.g. the Resonating Group Model Q, Antisym- 
metrized Molecular Dynamics [1C| . or the Microscopic 
Multicluster Model jl4j . Moreover, many-body type of 
calculations have also been performed on 9 Be: projected 
Hartree-Fock [1] , Shell Model 3 , Quantum Monte Carlo 
[IH , and ab initio no-core shell model [l5| . All of them 
are able to reproduce the low-lying energy spectrum and 
electromagnetic properties in fair agreement with the ex- 
perimental data available at the moment, though in gen- 
eral theoretical models predict more states than are seen 
experimentally. 

Somewhat surprisingly, the three-body decays of the 
9 Be resonances have been barely studied [HI, (28|. The 
inverse process may proceed through the resonances but 
non-resonant contributions are also important. Before 
facing this more complicated process, it is advisable to 
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FIG. 2: (Color online) Scheme of the different Jacobi coordi- 
nates for 9 Be. 



FIG. 1: Scheme of the experimentally known levels of Be 
below 6 MeV of excitation energy. 

get a good understanding of the resonance decay of 9 Be 
into aan. The experimentally known 9 Be states are 
shown in Fig.Q]for excitation energies below 6 MeV where 
all other particle thresholds than aan are closed. All 
these levels, apart from | , have a fairly large width, 
which makes it difficult to determine their properties. 
Many experimental efforts have been addressed towards 
this | "-state [H [H, HI] . They all agree in the small 
percentage of the decay taking place via the 8 Be ground 
state. So far no agreement has been reached regarding its 
main decay path, via 8 Bc(2 + ), 5 Hc(p), or direct. Much 
less is known about the decays of other low-lying reso- 
nances of 9 Be. Both, the \ and | + , seem to prefer to 
decay through 8 Be(0 + ), although especially the results 
for i need to be better established. 

The purpose of the present article is to report on com- 
prehensive calculations of the three-body properties of 
low- lying states in 9 Be. We give a survey of the short- 
distance structure of the resonances, their dynamic evo- 
lution across intermediate distances which often is re- 
ferred to as decay mechanism, and eventually reaching 
the large-distance asymptotics which reveal the complete 
set of momentum distributions of the fragments after de- 
cay. The two-dimensional energy correlations showed in 
Dalitz plots can be directly compared to the experimental 
data. This is the only information relating measurements 
with initial short-distance structure and decay mecha- 
nism. Extrapolations backward from data, therefore nec- 
essarily must be model dependent. We attempt to pro- 
vide an interpretation which is as physically meaningful 
as possible. 



II. THEORETICAL INGREDIENTS 

The decay of 9 Be into two a-particles and one neutron 
is obviously a three-body problem in the final state where 
the particles are far from each other. Furthermore the 
dominant structure at small distances is also of cluster 
nature for these low-lying resonances. The Hamiltonian 



for this cluster structure is then 

H = t^ i -^+ty^ j -r k) + V 3b , (1) 

i— i—1 

where p t and are momentum and mass of particle 
i, Pt is the total momentum, Vi is interaction between 
particle j and k. Here {i,j, k} is a cyclic permutation of 
{1, 2, 3}, and is a three-body potential depending on 
all three particle coordinates. It is convenient to substi- 
tute the position coordinates by the Jacobi coordinates, 
x and y defined as: 

x = J (n -r 2 ) 

V m 

/M12.3 ( m 1 r 1 +m 2 r 2 \ 

y = \ — - [ r 3 - , (2) 

where m is an arbitrary mass scale chosen as the nucleon 
mass, and p\ 2 and fJ-12,3 are the reduced masses. The 
Hamiltonian becomes 

H 2 3 
H = - — {Vl + Vl)+Y J V i {x,y) + V zh . (3) 
2m y ^— ' 

i=l 

In the present case we have two possible choices for Jacobi 
coordinates (see Fig. [5]), leading to different sets of (x, y)- 
coordinates. We use hyperspherical coordinates where 
the six coordinates are {p,a,9 Xl 4> x fi v ,(j)y}. The #'s and 0's 
refer to the directions of x and y, while a = arctan(x/y) 
and p = [x 2 +y 2 ) 1 / 2 are related to their sizes. Actually p 
is the only length coordinate, which describes the average 
distance from the center of mass. 

Resonances are computed within the formalism of com- 
plex scaling of the hyperspherical coordinates. This is 
particularly simple since only one coordinate, the hyper- 
radius p, has to be scaled. The Hamiltonian is complex- 
rotated, i.e. 

H e {p)=H(pJ e ). (4) 

We use the adiabatic expansion method and solve the 
Faddcev equations stepwise, that is first the angular then 
the (hyper) radial part [29|. The angular part of the 
Hamiltonian is first solved keeping fixed the value of p, 
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(5) 

where n labels the adiabatic components. Tq is the angu- 
lar part of the kinetic energy operator [2Sj . This provides 
a complete set of angular wave- functions, $„jm, that are 
employed to expand the total wave-function \& JM ; 
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where the p-dependent expansion coefficients, f n {p), are 
the hyperradial wave functions obtained from the coupled 
set of hyperradial equations [29| : 
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where V$b is a three-body potential used for fine-tuning 
and the functions P nn > and Q nn i are given for instance in 
[29j . The eigenvalues A„ in Eq.© enter in ([7]) as a part 
of the effective adiabatic potentials: 
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Resonances are usually understood as states with com- 
plex energy E = En — iT 'r/2, where Er is the energy of 
the resonance and T^j is the width. If we define now 



k = J2p,\E\/h 2 , with p being the reduced mass of the 
system, we then have that the asymptotic form of the 
resonance wave function is given by 
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where with (3 = i arctan(^r-). The first term oscillates 
while the second one diverges. After the complex scaling 
transformation (p — > pe lS ) the radial asymptotic behav- 
ior becomes 



for these energies. This is not necessarily correct at short- 
distances where all 9 nucleons (and their intrinsic struc- 
ture) may contribute in different (cluster) configurations. 

We use the (complex scaled) three-body model at 
all distances because the decay properties only require 
the proper description of the emerging three particles. 
Therefore, the angular eigenfunctions and eigenvalues in 
Eq.([5j) are complex, as well as all the terms entering 
in the coupled set of radial equations ([7]). The miss- 
ing information, if any, beyond the three-body struc- 
ture, is the initial structure at small distances. This 
piece, acting as a boundary condition, is parametrized 
through a short-range three-body potential of the form 
V ib = Scxp(-p 2 /b 2 ). 

Different three-body resonances correspond in general 
also to different three-body structures. As a conse- 
quence, the missing information going beyond the two- 
body correlations is in principle resonance dependent. 
The strength (and possibly also the range) in the three- 
body force is therefore adjusted individually to give the 
correct position of each of the resonance energies. This 
adjustment implies that the potential is angular momen- 
tum dependent but this is already a property of the two- 
body potential. The corresponding Hamiltonian for the 
three-body problem still exists as a non-local operator, 
but this feature is already present due to the angular mo- 
mentum dependence of the two-body interactions. It is 
then clear that this phenomenological fine-tuning is not 
arising from the presence of a genuine three-body inter- 
action. 

The energy dependence is all-decisive for decay prop- 
erties as evident in the exponential dependence of proba- 
bility for tunneling through a barrier. On the other hand, 
the three-body potential is assumed to be completely 
structure independent, and therefore only marginally in- 
fluencing the partition between different structures at 
large distances. However, this is an assumption which 
may be violated through the dynamic evolution from in- 
accurate initial small-distance boundary conditions pro- 
vided by the three-body potential. 
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III. SHORT-DISTANCE STRUCTURE 



which implies that when 8 > (3 the wave function goes 
to zero exponentially and the resonance can then be ob- 
tained as an ordinary bound state. True bound states 
remain unchanged under the coordinate rotation. 

For our particular case of 9 Be, the two-body interac- 
tions, Vi, are chosen to reproduce the low-energy scat- 
tering properties of the two different pairs of particles in 
our three-body system. We use the Ali-Bodmer a — a po- 
tential [3(| supplemented by the Coulomb potential be- 
tween a-particles, and the a-neutron interaction is taken 
from [3l|. The 9 Be- resonances are of three-body charac- 
ter at large-distances, since no other channels are open 



The short-distance structures are crucial for the ener- 
gies whereas dominating configurations at large distances 
are decisive for the observable decay properties. The con- 
nection between these two regimes contain information 
about the decay mechanism which therefore only is an 
observable effect precisely to the extent reflected in the 
final distributions. In other words, sensible theoretical 
models are indispensable to interpret the experimental 
results. In this section we extract and discuss short- 
distance bulk properties, that is effective potentials, en- 
ergies and partial- wave structure. 
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FIG. 3: (Color online) Real parts of the seven lowest adi- 
abatic potentials as functions of the hyperradius p for the 
1/2+ , 5/2", 1/2", 5/2+, 3/2+ and 3/2" low- lying resonances 
in 9 Be. 



A. Adiabatic potentials and energies 

Each of the adiabatic potentials entering in Eq. ([7) cor- 
responds to a specific combination of quantum numbers, 
i.e. partial-wave angular momenta between the particles 
in the different Jacobi systems. Usually only rather few 
adiabatic potentials are needed to achieve convergence. 
We show in Fig. [3] the real part of these adiabatic po- 
tentials as defined in (jHJ plus the three-body potential 
individually fitted for each spin and parity in order to 
reproduce the experimental resonance energies. 

We did not include the non-adiabatic diagonal parts in 
the figure (Q nn in (J7J , P nn = 0) because they usually are 
insignificant and in a sense more related to the coupling 
between the different radial potentials. The only excep- 
tion is the deepest potential for the l/2 + resonance. In 
this case the Q33 term is included, since this term is re- 
sponsible for the potential barrier that permits to hold 
the resonance (see [32[ for details). The imaginary parts 
of the potentials are small, oscillate, and go to zero for 
large values of p. They are mostly related to the widths. 

In Fig. [3] we observe that at small distances the lowest 
potentials have a pronounced well followed by a poten- 
tial barrier that are responsible for the bound states and 
the resonances. We have two attractive potentials for 
angular momentum and parity J T = 3/2"; the deepest 
one supports the ground state (bound) while the other 
one supports the higher resonance. For other J T values 
only the lowest potential exhibits an attractive region at 
small distance. All the other potentials are repulsive at 
all distances. 

At large distances, the lowest potential is in all 
the cases approaching the 8 Bc(0 + ) resonance energy of 
~ 0.1 McV. Its angular structure corresponds to the two 



TABLE I: Calculated and measured energies Er (in MeV) 
and widths Tr (in MeV) of the 9 Be resonances for different 
J*. The experimental values (labeled "exp") are from [33|] 
and the calculated results (labeled "th") are obtained with 
the three-body interaction parameter S (in MeV) (the range 
is taken b = 5 fm in all the cases). The energies are measured 
from the ctan threshold. 
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alpha particles populating the 8 Be(0+) resonance, while 
the remaining neutron is far away and described through 
the radial equation. These specific potentials are labeled 
as number n = 3 in 1/2+, n — 7 in 5/2", n = 1 in 1/2", 
5/2+ and 3/2+, and n — 2 in 3/2". They characterize 
then a decay mechanism where the neutron first is emit- 
ted and this two-body a — a resonance is populated and 
subsequently decaying. 

The complex scaling of the hyperradius leads to a 
Hamiltonian with complex solutions vanishing exponen- 
tially at large distances precisely as ordinary bound state 
wave functions. The real and imaginary parts of the com- 
plex three-body energy are, respectively, the resonance 
energy Er, and — Tr/2, where Tr is the width of the 
resonance. The computed results are collected in table U 
together with the known experimental values [33j. The 
three-body strength is adjusted to give the correct en- 
ergy position for all J* except for 3/2", where we choose 
the same values for both bound state and resonance as 
in ref. 0. 

The energies are given relative to the a — a — n break- 
down threshold, at 1.57 MeV above the ground state. 
We have found one bound state, which corresponds to 
the 9 Be ground state, and six resonances below 6 MeV 
of excitation energy. The ground state has J 71 " = 3/2" 
and the resonances l/2 ± , 3/2+ and 5/2+. Those states 
are most likely to contribute to processes bridging the 
A = 5, 8 instability gaps in nuclear synthesis in suitable 
astrophysical environments (35| . We keep the range of the 
three-body potential at b — 5.0 fm for all J v , while ad- 
justing the strength to place the resonance energies (and 
bound state) at the desired measured position. Thus, we 
did not attempt to reproduce the widths. 

Several features are interesting in table |TJ First 
the three-body potentials have very moderate strengths 
which were found to reproduce the real part of the mea- 
sured energies. This is fine-tuning and indicates strongly 
that the dominating structures in fact really are three- 
body clusters. It is then significant that all widths, ex- 
cept for l/2 ± and the very narrow 5/2" state, are larger 
than the corresponding measured values. This is consis- 
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tent with a many-body configuration at small distances 
which would decrease the branching ratio of decay into 
the investigated three-body cluster structure. A possible 
quantification of this deviation is in terms of preforma- 
tion factors expressing that only part of the complete 
wave-function describes the three-body cluster. Again 
the deviation amounts to about factors of two in agree- 
ment with only smaller contributions from a many-body 
structure. A smaller range compensated by a slightly 
larger strength to leave the energy untouched would de- 
crease the width towards the measured values. 

However, the 3/2 _ state is another exception in ta- 
ble |U The computed value of width is smaller than the 
experimental table value. The computed energy is also 
below measurements by about 1.4 MeV. An attempt to 
increase the computed energy by this amount with use 
of a repulsive three-body potential of range 5 fm immed- 
uately cause the resonance to disaapear into the contin- 
num corresponding to an energy above the barrier. This 
shows that the width in the model very quickly becomes 
very large, and exceeding the experimental table value. 
Either the model is missing an important ingredient for 
this state or its width should be substantially larger. 

The l/2 ± states are an apparent exception in table U 
where the experimental widths are larger than the cal- 
culated values. For l/2 + this discrepancy has caused a 
good deal of trouble. In a recent investigation [32| this 
conundrum was explained as a genuine three-body effect 
where the resonance structure changes from dominantly 
5 He plus an a-particle at small distances to 8 Be plus a 
neutron at large distances. The large measured width is 
in fact obtained from an assumption of two-body char- 
acter and sequential decay in the i?-matrix analysis of 
the photo dissociation cross section [23| . It is remarkable 
that a much smaller width consistent with [32j was ob- 
tained already many years ago H . The large experimen- 
tal width of 1/2" should be reevaluated since we suspect 
that the i?-matrix parametrization and the sequential de- 
cay channel is used too strongly in the extraction from 
the data analysis. This was argued for 1/2+ in (32[. 



B. Partial waves 

The different two-body components of the three-body 
system are constrained by the total angular momentum 
and parity of each state. For example, in the first Ja- 
cobi (see Fig. [2]) the two a-particles must couple to an 
even orbital angular momentum £ x . A neutron with even 
(odd) angular momentum, £ y , will give a positive-parity 
(negative-parity) state. The orbital angular momentum 
£ y couples to the spin of the neutron to the angular mo- 
mentum j y . In the second Jacobi set £ x can be either 
even or odd, and couples to the spin of the neutron to 
the angular momentum j x . 

We choose our partial-wave components taking into 
account these selection rules. In the present case conver- 
gence is achieved with a number of partial waves between 



TABLE II: Components included for each state of 9 Be 
relative to the first Jacobi set (see Fig. [5J . q labels the set of 
quantum numbers, £ x , £ y and j y (coupling between £ y and the 
spin of the neutron) are the angular momenta relative to x 
and y Jacobi coordinates. K max is the maximum value of the 
hypermomentum and W gives the probability in % for finding 
these components in the resonance. Only the components 
contributing more than 1 % are shown. 
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10 and 30, depending on the resonance. The accuracy is 
optimized by choosing a large value for the hypermomen- 
tum K for the large contributions. Unfortunately, the 
higher the value of K maxi the larger becomes the total 
number of basis states. Therefore, K max must be chosen 
carefully for each partial wave, trying to achieve accuracy 
while keeping the number of basis elements as small as 
possible. 

Tables [TT] and IIIII show, for the first and second Ja- 
cobi sets, the contributions W to the total J* wave func- 
tions from those components contributing more than 1%. 
These contributions are well defined for the complex ro- 
tated resonance wave function, since it behaves asymp- 
totically like a bound state. The maximum value of the 
hypermomentum K is also given for each component. 
The computed values of W are given in the last column 
of the tables. The wave functions are located at relatively 
small distances, and the contribution from the different 
components (obtained after integration of the square of 
the wave function over all the hyperangular variables) 
contain therefore information mainly about these bulk 
structures. The decay properties are contained in the 
large-distance tails, whose partial wave content can be 
entirely different, as discussed in details in the next sec- 
tion. 

The lowest resonance, J 71 " = i + , is located only 18 keV 
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TABLE III: The same as table[TT]for the second Jacobi set (see 
Fig-©. The angular momentum j x results from the coupling 
between £ x and the spin of the neutron. 
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above the two-body 8 Be narrow ground state resonance 
at 918 keV. In the first Jacobi coordinates this state is 
entirely described as s-waves between the a-particles and 
therefore also between their center of mass and the neu- 
tron. The interesting structure is seen in the two other 
identical Jacobi coordinates where the structure changes 
abruptly from a— neutron p 3 / 2 to Si/ 2 configurations at 
around 10 fm, see [32[. The bulk part of the resonance 
structure found at small distances then roughly amounts 
to equal parts in each of these partial waves. 

The next resonance with J 71 = | is very narrow due 
to the large barrier in the dominating partial wave of 
(£x,£y) = (2, 1) in the first Jacobi and (1,2) in the sec- 
ond set of Jacobi coordinates. This can be described as 
8 Bc(2 + ) or 5 He(p 3 / 2 ), respectively, but it is in fact the 
same state in different coordinate systems. It is therefore 
not meaningful to distinguish between these configura- 
tions unless also spatial distributions are included in the 
distinction (28j . 

The ^ resonance is a result of the 5 He p-wave attrac- 
tion combined with orbital angular momentum coupling 
to 2 of the last a-particle. Only the corresponding adi- 
abatic potential is really attractive. This configuration 



translates to (£ x ,£ y ) = (0, 1), (2, 1) in the first Jacobi co- 
ordinates where only even £ x are allowed. 

The | + resonance is dominated by a combination of 
5 He(p3/ 2 ) and 8 Be(0 + ,2 + ). Only one of the adiabatic 
potentials is really attractive and in fact not very deep. 
This state is important at moderate temperatures for 
photo dissociation and three-bod y r ecombination from 

the continuum via £T-transitions |35| . 

3 + 

The next resonance, | , is higher. Its structure is 

similar to the | + state in the first Jacobi where (£ x ,£ y ) = 
(0,2), (2,0) are roughly interchanged in the two states. 
In the second Jacobi system the 5 He(p 3 / 2 ) structure also 
has a relevant, although not dominant, contribution. 

The last resonance, | , has (£ x ,£ y ) = (2, 1) and (1,2) 
in the first and second Jacobi sets, respectively. This 
is reflecting a combination of the influence of the inter- 
actions related to the 5 He(p 3 / 2 ) and 8 Be(2 + ) two-body 

resonances. The similarity to the | state is striking, ex- 
cept for the larger width arising from a higher excitation 
energy. 



IV. LONG-DISTANCE STRUCTURE 

Resonances may be populated at small distances via 
beta-decay or some specific reactions, but the products 
after the resonance decay reflect the behavior at large 
distances. The short and large-distance structures are 
related through the quantum mechanical solution, and 
the configurations sometimes change dramatically with 
the hypcrradius. This connection from small to large 
distances is therefore crucial for the interpretation of the 
decay mechanism and the measured results. We shall first 
show the dynamic evolution of each resonance configura- 
tion, and afterwards show the momentum distributions 
of the fragments as Dalitz plots with the full information. 



A. Dynamic evolution 

As mentioned in section UlI A[ at large distances, the 
lowest adiabatic potential is, for all the resonances, ap- 
proaching the 8 Bc(0 + ) resonance energy of ~ 0.1 MeV. 
Its angular structure corresponds to the two alpha parti- 
cles populating the 8 Be(0 + ) resonance, while the remain- 
ing neutron is far away. In other words, for large values 
of p, the configuration of this potential in the first Jacobi 
set approaches £ x = 0, j y = J, and £ y has to be one of 
the J ± 1/2 values in order to produce the correct parity. 

If this state is populated at large distances, where all 
couplings to other adiabatic potentials have vanished, the 
decay can be described as sequential via the 8 Be two- 
body ground state. Such a decay has a special role 
because it is favored by a very low energy with non- 
vanishing coupling to other potentials for all the reso- 
nances. Figs. 2] and [5] show the partial- wave decomposi- 
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25 50 75 25 50 75 25 50 75 
P(fm) 

FIG. 4: (Color online) Contribution, as a function of p, of the 
different partial waves in the first Jacobi set (labeled by q as 
in table |TT| to the adiabatic eigenfunction $ 9 (p, fi) related to 
the ground-state structure of 8 Be at large distance. 




25 50 75 25 50 75 25 50 75 
P(fm) 

FIG. 5: (Color online) The same as in Fig. [4] but for the 
second Jacobi set. q labels the set of quantum numbers as in 
table HTI1 



tion, as a function of the hypcrradius, for this adiabatic 
component. 

Not surprisingly, the dominant partial wave in the first 
set of Jacobi set is the 8 Be(0 + ) structure for p values 
beyond about 10 fm, see Fig. 2] The same simple struc- 
ture does not appear in the second Jacobi coordinates 
as seen in Fig. [51 At short-distances, around 20 fm, one 
of the components gives most of the contribution, but 
this structure is not maintained at large-distances where 
we observe a very fragmented partial-wave decomposi- 
tion. The reason is of course that the transformation of 



the l x = state in the first Jacobi set into the second 
one results in contributions from many different angular 
momentum components. 




25 50 75 25 50 75 25 50 75 
P(fm) 



FIG. 6: (Color online) Contribution, as a function of p, of the 
different partial waves in the first Jacobi set (labeled by q as 
in table [TTJ) to the dominant adiabatic eigenfunction <& q (p, Q) 
different from the one related to the ground state structure of 
8 Be at large distance. 




25 50 75 25 50 75 25 50 75 
p(fm) 

FIG. 7: (Color online) The same as in Fig. [6] but for the 
second Jacobi set. q labels the set of quantum numbers in 
each component as in table [TTT1 

The sequential decay via the 8 Bc ground state leads 
to very simple momentum distributions derived from the 
two-body character where energy and momentum con- 
servation fully determine the final state. The remain- 
ing part of the decay proceeds through other adiabatic 
components. We show in Figs. [5] and [7] the partial- wave 
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decomposition as function of hypcrradius for the most 
contributing of these other components. 

Each of the J T states presents its own features. In all 
cases the variation from small to large distances is sub- 
stantial and sometimes dramatic. The structures always 
converge at large distances. In the first Jacobi (Fig. [5]), 
one of the partial waves absorbs after convergence al- 
most all the contribution for the resonance of 1/2+, l/2~, 
5/2+ and 3/2 + . The corresponding partial waves have 
£ x = 0,0,2 and 2, respectively. They are, therefore, re- 
lated to + and 2 + states in 8 Bc. In the 5/2" resonance 
there are two partial waves contributing significantly, but 
one of them is dominating with i x — 2 and thus related 
to a 8 Be(2 + ) structure [28|]. In contrast, the 3/2~ reso- 
nance has two components with l x = 2 that are equally 
important at large values of p. 

We show in Fig. [7] the partial- wave decomposition for 
these resonances in the second Jacobi system. The 1/2+ 
state presents a mixture of three components contribut- 
ing significantly at large distances. This is somewhat 
analogous to the transform of the 8 Be ground state with 
many partial waves at large distances. Here we only find 
three which therefore also emphasizes that even with the 
same quantum numbers (all s-waves) the structure can 
be very different. The 5/2 ~ and l/2~ states reveal struc- 
tures where each converges to quantum numbers identical 
to those of 5 He(p3/2) and 5 Hc(p 1 / 2 ) respectively. Again 
this does not imply that these are the decay channels, 
only the quantum numbers are the same. In 5/2 + and 
3/2 + the dominant components at large distance have 
I x = as for the 8 Be ground state but the behavior dif- 
fers very much from those of Fig. [5] The last state of 
3/2~ presents two roughly equal contributions both with 
(£x,£y) = (1, 2) very much like in the first Jacobi system. 
It is like 8 Be(0 + ) is replaced by 5 He(p-state). 



B. Momentum distributions 

The decay mechanisms depend on the resonance prop- 
erties and they are conventionally called either sequential 
via a given two-body structure or direct decay to the con- 
tinuum or perhaps a mixture of these possibilities. The 
process is sequential when the measured kinematics re- 
veals that one particle is first emitted and subsequently 
the remaining structure decays into two fragments in- 
dependent of the emission of the first particle. Com- 
binations of such decay channels form the basis for the 
R- matrix analyses of experimental data (36j | . This formu- 
lation becomes dubious when the intermediate two-body 
structure falls apart on the same time-scale as the first 
emission. The process is then better described as a gen- 
uine three-body decay. This does not prevent analyses 
in terms of several two-body decay channels. The two 
different formulations may still be completely identical 
provided the two different sets of basis functions span 
the same space. The two formulations merely differ in 
the choice of basis [37| ■ 



The present case is in one way rather simple since all 
the resonances can decay sequentially via 8 Be(0 + ), which 
is a long-lived stable structure surviving long time after 
emission of the neutron. Thus the sequential decay mech- 
anism through this state is not controversial and easily 
separated kincmatically in experiments. One of the adi- 
abatic components is related to the 8 Be+n structure and 
approaches the energy of the 8 Bc(0 + ) resonance. This 
component describes the sequential decay contribution 
through this channel. Extension of this picture to se- 
quential decays through 8 Be(2 + ) or 5 He(p 3 / 2 ) is an in- 
vitation for difficulties, since these channels are broad 
(short-lived) resonance structures, not easily separated 
from the background continuum, and furthermore not 
even orthogonal contributions. This is more reasonably 
described as a direct decay. 

The technique involved is described in [38| where the 
large-distance asymptotic behavior of the radial wave 
functions are shown to give the "branching ratio" for such 
sequential decay. We have calculated this fraction of de- 
cay for each of the 9 Be resonances. The result is given 
in table IIV1 The decays of \ are found to be predom- 
inantly sequential. In | + both mechanisms are compa- 
rable whereas the direct decays dominate for the other 
three resonances. The comparison to measured branch- 
ing ratios is rather favorable in view of the uncertainties 
for broad resonances and the different methods of extrac- 
tion. 

The uncertainties arc especially emphasized by con- 
sidering the g state which often is quoted as predomi- 
nantly decaying through the 8 Be ground state 0,13 in 
agreement with our result. This is intuitively appealing 
since the alternative channels of 8 Be(2 + ) and 5 He(pi/ 2 ) 
are rather high-lying. More recently also contributions 
through such channels are extracted from experimental 
analysis [25|, [27| although given with reservations and un- 
certainties. Furthermore, the beta-feeding, the width, 
and the decay channel are linked together for broad res- 
onances in data analysis [2l[. We conjecture that the 
width should be smaller than the measured value in ta- 
ble U and the predominant decay channel is through the 
ground state of 8 Be. 

For the sequential channels, the resulting momentum 
distributions are easily found, since the first emission im- 
mediately provides the energy of the particle in the three- 
body center of mass system. The following decay is again 
given by one energy in the center of mass system of the 
remaining two particles. 

The momentum distributions for direct decays into 
the continuum can now be found by excluding the se- 
quential contribution, that is the part of the wave func- 
tion residing in the 8 Be ground state at large distance. 
Again we have to calculate, as accurately as possible, 
the large-distance asymptotics of the wave function. The 
technique, described in [4^, EH , is based on finding the 
Zeldovic regularized Fourier transform of the coordinate 
state wave resonance function. The result is directly 
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TABLE IV: 9 Be resonance excitation energies, energies above the aan threshold, and, for each resonance, estimated amount 
of computed and observed sequential decay via 8 Be(0 + ). 

r E aan (MeV) E re3 (MeV) Theo. (%) Exp.(%) [27] Exp. (%)[39, 40] Exp. (%) [41] 



5 + 



1.68 


0.11 


100 






2.43 


0.86 


3 


6±1 


7± 1 


2.82 


1.25 


90 


32 ± 15 


100 


3.03 


1.46 


53 


46 ±20 


87 ± 13 


4.69 


3.12 


1 


16 ±2 




4.22 


2.65 


29 







100 



comparable to measured distributions. It is worth em- 
phasizing again that the only link from the asymptotic, 
measurable distribution, to the small-distance structure 
is via theoretical models [28| . 

We compute the distributions by Monte Carlo simu- 
lation. We first generate randomly a large number of 
events, each of them consisting of three four-momenta 
relative to our three decaying fragments. The sum of 
their center-of-mass energies must equal the resonance 
energy. The weight of each set of momenta is the 
absolute-squared wave-function at large distance. The 
resulting energy distributions are shown in Figs. |S] and 
[9] for a-particles and neutrons, respectively. We give the 
energies in units of their maximum values for each case, 
i.e. 5/9E res for the a's and 8/9E res for the neutrons. 

The distributions all necessarily have peaks, since they 
start with zero and return again to zero at maximum en- 
ergy. However, they can have more than one peak, and 
each of them has an individual position and width. The 
5/2~ and 3/2~ resonances are smooth with one peak for 
both neutrons and a-particles. In both cases, the neutron 
energies peak below and the a-particle above half of their 
respective maximum values. This means a tendency to 
emit a-particles in essentially opposite directions while 
leaving the remaining neutron in the middle with rela- 
tively little energy. This is only a tendency and the full 
distributions require detailed computations. Still it is 
indicative for this part of the process. 

For the l/2~ decay the a-particle show up with a broad 
distribution on the low-energy side whereas the neutron 
appears on the high-energy side. This resembles the se- 
quential decay through the 8 Be ground state where a- 
particles end up with only little energy. However, the 
present decay has to proceed through an orthogonal adi- 
abatic potential which reveals itself by the low-energy 
node in the distribution of the neutron energy. 

The 3/2 + and 5/2 + resonances both produce neutrons 
and a's with tendencies to be respectively on high and 
low-energy sides like for the 5/2~ and 3/2~ resonances. 
However, in the positive parity cases additional peaks 
appear in both distributions, again a signal of an excited 
state. These cases are otherwise not very similar and the 
distributions are very broad each extending across from 
high to low-energy side and vice versa. 

These distributions can be suggestive and deceiving. 
The momenta are distributed among all the three parti- 



cles which is the reason for the continuous distributions 
in the first place. However, this also means that a kine- 
matically complete description for a given conserved total 
resonance energy requires energies of two particles at the 
same time. This information is contained in the two- 
dimensional energy correlations known as Dalitz plots 
which were introduced by R.H. Dalitz in 1953 to study 
decays of K- mesons [3] ■ These correlation diagrams pro- 
vide an excellent tool for studying the dynamics of three- 
body decays. The technique has recently been picked up 
and applied in studies of nuclear fragmentation processes 
[45j . In simple two-body decays the angular distribution 
of the emitted particles carries the signature of decay- 
ing angular momentum and parity. The Dalitz plots are 
generalizations to three-body decays and it is natural to 
use the plots in attempts of experimentally assigning spin 
and parity to the decaying resonances, [371 ]. 

To establish the connection to measured distributions 
we computed Dalitz plots for a-particles and neutrons af- 
ter the decay of 9 Be resonances. We use the same Monte 
Carlo technique as for the individual particles energy dis- 
tributions. To facilitate comparison with the experimen- 
tal results from (27j we plot the a — n relative energies 
on the x and j/-axes, i.e. 

t-^ |Pa — n\ L^an ^an 

E a - n = — , p a _ n = —P a ~ —Pn , (11) 

where p Q _ n is the relative momentum. The results are 
shown in Fig. [TO] We first observe that all the distribu- 
tions are symmetric with respect to interchange of the 
axes. This is necessary and reflects that the wave func- 
tions are symmetric for the identical bosonic a-particles. 

The graphs corresponding to | and | , are very sim- 
ilar to each other. None of them exhibits any points or 
regions of zero probability, except on the confining en- 
velope defined by energy conservation. This means that 
symmetry, angular momentum and parity of these struc- 
tures, (£ x ,£ v ) = (1,2), (2,1) (tables HI and |in]), allow 
emission in all directions and with all energy partitions. 

The probability increases towards higher a energies, 
which corresponds to smaller relative energies since the 
neutron is much lighter than the a-particle. This is due 
to the Coulomb repulsion and the tendency to choose a 
decay path where the neutron is left in the middle as ob- 
served in the one-dimensional energy distributions, sec 
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FIG. 8: (Color online) Q-particle energy distributions for a) the 5/2 -resonance of 9 Be at 2.43 MeV of excitation energy (or 

0. 86 MeV above the aan threshold), b) the l/2"-resonance of 9 Be at 2.82 MeV of excitation energy (or 1.25 MeV above the 
aan threshold), c) the 5/2 + -resonance of 9 Be at 3.03 MeV of excitation energy (or 1.46 MeV above the aan threshold), d) the 
3/2 + -resonance of 9 Be at 4.69 MeV of excitation energy (or 3.12 MeV above the aan threshold), e) the 3/2 "-resonance of 9 Be 
at 4.22 MeV of excitation energy (or 2.65 MeV above the aan threshold). The energies are divided by the maximum possible, 

1. e. 5/9E reB . The sequential decay via 8 Be(0 + ) has been removed in all the cases. 




0.5 1 0.5 1 0.5 1 0.5 1 0.5 1 

E|-/En,max E n /E n max E n /E nimax E n /E n max E n /E n max 



FIG. 9: (Color online) The same as Fig. [8] for the neutron energy distributions. The energies are divided by the maximum 
possible, i.e. 8/9E res . 



Figs. [5] and [5] Distributions for both states compare well 
with the experimental plots from J27|. Moreover, the | 
distribution is very similar to the measured ones pub- 
lished in [26|, and investigated theoretically in detail in 

The distributions for | , | + , | + exhibit much more 
structure and all have zero probability regions as reflected 
in the nodes or minima of the one-dimensional distribu- 
tions. For i we find a striking similarity with the mea- 
sured distribution in p7j at the excitation energy window 
at 2.8 MeV. The very low probability bands at lower and 
higher relative energies are found both places. The dif- 
ferent projection on the neutron energy axis in Figs. [9] 
resulted in a node at small neutron energy, presumably 
corresponding to a cut along low energy small probability 
region in Fig. ITUl 

For both | and | the computed distributions have 
lots of structure whereas the measurements show more 
smooth distributions without much resemblance to cal- 
culations. The explanation for these discrepancies is still 
to be found. 



V. SUMMARY AND CONCLUSIONS 

The hyperspherical adiabatic expansion method, com- 
bined with complex scaling, is used to compute the ener- 
gies and widths of 9 Bc low-lying resonances. We describe 
them as three-cluster resonances (aan). Realistic short- 
range nuclear interactions as well as Coulomb interac- 
tions are included in the computations. To reach high 
accuracy we use a large hyperharmonic basis for each 
angular eigenfunction, accurate large-distances, outgo- 
ing waves of radial wave functions and, if possible, the 
correct energy of the three-body resonance obtained by 
tuning the three-body potential. 

We find one bound state (| ) and six resonances below 
6 MeV of excitation energy in agreement with experimen- 
tal information. Spins and parities of the resonances are 
i , | and | . The small-distance properties of the 
adiabatic potentials determine energies, while barriers at 
intermediate distances are crucial for the widths, and the 
large-distance structure of the resonances are decisive for 
the momentum partition between the three particles in 
the final state after decay. 

The structure of the resonances are obtained as differ- 
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FIG. 10: (Color online) Dalitz plots for a) the 5/2~-resonance of 9 Be at 2.43 MeV of excitation energy (or 0.86 MeV above 
the aan threshold), b) the l/2~-resonance of 9 Be at 2.82 MeV of excitation energy (or 1.25 MeV above the aan threshold), 
c) the 5/2 + -resonance of 9 Be at 3.03 MeV of excitation energy (or 1.46 MeV above the aan threshold), d) the 3/2 + -resonance 
of 9 Be at 4.69 MeV of excitation energy (or 3.12 MeV above the aan threshold), e) the 3/2~-resonance of 9 Be at 5.59 MeV of 
excitation energy (or 4.02 MeV above the aan threshold). We plot on the axis the a-n relative energy in MeV. The sequential 
decay via 8 Be(0 + ) has been removed in all the cases. 



ent combinations of angular momenta of the two-body 
subsystems. The configurations are determined by the 
interactions leading to observed low-lying resonances of 
the subsystems, i.e. + ,2 + for 8 Bc and p 3 / 2 P1/2 f° r 
5 He. The detailed configurations of the three-body res- 
onances are extracted, their energies fine-tuned via the 
three-body potential, and their widths computed. 

We compute the possibly substantial dynamic evo- 
lution of the resonances as functions of hypcrradius. 
The large-distance asymptotic structures are via Fourier 
transformation directly related to the momentum distri- 
butions of the fragments after the three-body decay. We 
determine the fraction decaying via the ground state of 
8 Be in a sequential decay. The agreement with measure- 
ments is rather good in view of the uncertainties related 
to broad resonances and different theoretical and exper- 
imental definitions and methods. 

The remaining part is described as direct decay to 
the three-body continuum. We present the computed 



momentum distributions of neutrons and a-particles for 
each of the resonances. These observable distributions 
are results of the dynamic evolution, and open to exper- 
imental tests. We compare with the available data, and 
find remarkable similarities except for the 5/2 + , 3/2 + res- 
onances where the theory gives much more structure than 
found in the energy windows selected in the experiments. 
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